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ABSTRACT 

We present a new measurement, 0.85- 3. 5Zq, of the metallicity of high velocity cloud (HVC) Complex 
M by analyzing ultraviolet spectroscopic observations of the blazar Mrk 421 taken with the Cosmic Origins 
Spectrograph on the Hubble Space Telescope and the Far Ultraviolet Spectroscopic Explorer. Although an 
HVC at Vlsr = -131 km s" 1 is not visible in 21 cm emission (logA/ni < 18.38; 3a), it is detected in ultraviolet 
absorption lines of C n, N i, O i, O vi, Si n, Si m, Si iv, Fe n, and H i. By referencing velocities to the intermediate 
velocity cloud at -60 km s" 1 and jointly analyzing H i absorption from high-order H i Lyman lines, we measure 
logA/m = 16.84^| (la) in the HVC. Comparing H i and O i, we find an HVC metallicity [0/H]= 0.32^;g. 
Because the sight line passes ~ 4° from the HVCs in Complex M, the detected HVC may represent the highest 
velocity component of the Complex, and our measurements provide a lower limit to its metallicity. The high, 
possibly super-solar metallicity, together with the low distance, z < 3.5 kpc, above the Galactic plane suggest 
that Complex M is condensed returning gas from a Galactic fountain. 

Subject headings: ISM: clouds — ISM: abundances — ISM: structure — ultraviolet: ISM 



1. INTRODUCTION 

The gaseous clouds surrounding the Milky Way are ob- 
served at near-zero velocity, as well as in intermediate ve- 
locity clouds (IVCs) and high velocity clouds (HVCs). The 
latter are defined as interstellar gas with absolute veloci- 
ties relative to the local standard of rest (LSR) between 
|Vlsr| = 30-90 km s" 1 and |V LS r| > 90 km s" 1 , respec- 
tively. These clouds are inconsist ent with the si mple model 
of Galactic rotation (see review bv lWakkerll200Tl) . The IVCs 
and HVCs are detected not only in neutral and mild ly ion- 
ized species like H i, C n, O i, S n, and Si n ( Murph 7et al.l 
[1991 IWakker & van Woerder3[l997l: ICollins et alJl2003lT but 
also in highly ionized species like C iv, O vi, N v, a nd Si iv 
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(Sembach et al. 
iLehner & Howk \ 

The IVCs have near-solar metal abundances, 



ICollins et all 12001 IFox et al.l [2006; 

I I2010I) . 

are found 

within several kpc of the Galactic plane (e .g., IWakker et al.l 
1996: iRvans ellffl[l997l: iRkhter et al.l l2001). and are consis- 
tent with the Galactic fountain interpretation. However, the 
origin of HVCs is still poorly understood. The proposed 
scena ri os include the Ga l actic fountain dShapiro & Field! 
119761; iHouck & Bregmanl 119901) . co oling Galactic halo 
gas falling onto the Galactic plan e (IWakker et al.l 119991; 
ICollins et all 120031; IShull et al.1 120091) . warm-hot intergalac- 
tic medium (WHIM) in the Local Group (Bli tz et ail 
119991; iBraun & Burtonll 19991) . and combinations of all three 
(Tripp et al. 20031). These competing models predict distinct 
scale sizes and metallicities of HVCs, ranging from several 
kpc with near-solar abundances (fountain model), to tens of 
kpc with subsolar abundances (halo model), to hundreds of 
kpc with primordial abundances (the WHIM model). Clearly, 
distances and metallicities are key measurements for con- 
straining the origin of HVCs. 
The HVC Complex M lies at the highest Galactic latitude 
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FIG. 1. — 21 cm emission observations toward Complex M region. Grey- 
scale image shows HVC in range Vlsr = - 80 to -130 km s -1 . Black contours 
indicate brightness levels 0.2, 0.5, and 1 K, and red contours indicate IVCs 
with brightness levels 0.5, 1, and 2 K. Three sight lines of interest and core 
regions of MI, Mil, Mill an d their co r respon ding IVCs are also marked. Fig- 
ure is similar to Figure 5 of Wakker (2001) but using observations of LAB 
H I survey. 

among known HVC complexes, covering b = 45° to 70° and 
longitudes from i = 130° to 200°, over a velocity range of 
-125 < vlsr < -85 kms" 1 (Hulsboschj[l968|). It is com- 
posed of several distinct HVC clouds (MI, Mil, and Mill) 
that are apparently superposed with several intermediate ve- 
locity clouds in the sky, including the IV Arch (IV2, IV6, and 
IV3: Fig. [Tl iGiovanelli etall [19731; I Wakker & van Woerdenl 
[19911 iKuntz & DanivM T996). In a comparison of spectro- 
scopic observations of two Galactic halo stars HD 93521 and 
BD+38°2182 (Fig. CD) observed w i th the International Ultra- 
violet Explorer (IUE). iDanly et alJ (119931) found high velocity 
wings at Vlsr < _ 85 km s" 1 in O i, Si m, and C n absorption 
in spectra of BD+38°2182 (z « 4.4 kpc) but not in HD 93521 
(z ~ 1.5 kpc). They estimated the distance of Complex M 
(Mil and Mill) to lie between 1.5 < z < 4.4 kpc above the 
Galactic plane. Rv ans et al] 0997) made further differential 
optical observations of stars in the region and revised the dis- 
tance (z = 3.5 kpc) to BD+38°2182. They argued that the ab- 
sence of high-velocity absorption along the HD 93521 sight 
line could be due to small-scale spatial variation of the Com- 
plex, and suggested decreasing the lower limit to its distance. 

The metallicity of Complex M is much less certain. Com- 
bining a much deeper 21 -cm emission limit, log A/hi < 17.43 



2 



Yao et al. 



from the Jodrell Bank 76m, with O i absorption measure- 
ment s by iDanlv et alj ([1993b along the BD+38°2182 sight 
line, iRvans et al.l d 1997 ) est imated an oxyg en abundance > 
0.5Z Q for the Mill cloud. IWakkeri (120011) constrained the 
metallicity of the MI cloud to be 0.4-1.8 Z Q based on S n 
and H i emission data in the core region of MI; the large range 
arises from uncertainty in the ionization correction. Because 
of its short distance and possible high metallicity, Complex 
M is believed to be part of the IV Arch dKuntz & Danlvll 19961; 
lWakkeril200l . 

In this Letter, we analyze an HVC along the sight line to 
the blazar Mrk 421 and set a tigh t lowe r limit to the metallic- 
ity of Complex M. ISavage et all d2005h detected the HVC in 
H i, C n, and C m absorption lines in the FUSE spectrum, but 
they focused on the origin of the extended red wing of O vi ab- 
sorption. Throughout, we quote errors at 1 a and upper limits 
at 3 a confidence levels. We adopt atomic data fr om Morton 
(2003) and solar abundances from lAsplund et all (120091) . in 
particular, log(X/#) + 12 = 8.43, 7.83, 8.69, 7.51, 7.12, and 
7.50 for C, N, O, Si, S, and Fe, respectively. 

2. OBSERVATIONS AND DATA REDUCTION 

Our data come from far ultraviolet spectroscopic ob- 
servations of the background target Mrk 421 (l,b = 
179.83°, 65.03°; z = 0.03 ) obtained with the C osmic Ori- 
gins Spectrograph (COS; Osterma n et al.l BOllb aboard the 
Hubble Space Telescope (HST) and with the Far Ultravio- 
let Spectroscopic Explorer (FUSE). Mrk 421 is one of the 
Guaranteed Time Observation (GTO; PI: Green) targets of 
COS, which has two medium-resolution gratings G130M and 
G160M, covering wavelengths 1150- 1450 A and 1405-1775 
A with a spectral resolution ~ 18,000. The FUSE data cover 
905- 1187 A with a resolution - 20,000. COS observed 
Mrk 421 on 2009 November 24 (PID 1 1520) with total expo- 
sures of 1.7 ks (G130M) and 1.8 ks (G160M). FUSE observed 
the source on 2000 December 1 (PID P101) and 2003 January 
19 (PID Z010) with a total exposure of 85.0 ks. 

The observations and calibrated data were retrieved from 
the Multimission Archive (MASlQ) at the Space Telescope 
Science Institute. The COS data were further processed 
following the procedure described in Danforth e t al.l (12010b . 
Only the FUSE LiFl a and SiC2a da t a were used, reduced fur- 
ther as described in Danforth et al. (2006). Individual expo- 
sures were weighted by their exposures and coadded to form 
final spectra, whose signal-to-noise ratios per resolution el- 
ement range from 25-35 (COS) and 5-15 (FUSE). We find 
< 5% flux fluctuation in the COS spectra arising from unre- 
moved small fixed-pattern features. These fluctuations, com- 
pared to systematic uncertainties in wavelength calibration, 
are not important in measurements described below. 

The wavelengths of these spectra were calibrated with 
respect to the H i 21 cm emission data from Lei- 
den/Argentine/Bonn (LAB), which have an angular resolu- 
tion of ~ 0.6°, spectral resolution of 1.3 k m s" 1 , and rms 
brightness-temperature noise of 70-90 mK (Kalberla et al. 
2005). Four LAB pencil beams adjacent to the Mrk 421 
sight line were extracted, averaged, and weighted by the in- 
verse square of their angular separations from the sight line. 
An IVC component is clearly resolved from the zero-velocity 
component in the final emission spectrum. We use a Gaussian 
fit to the IVC and obtain its LSR velocity, column density, and 



TABLE 1 
HVC Measurement 



ions 


b 


log/V 




(kms -1 ) 


(N in cm -2 ) 


Si II 


13.9 ±1.1 


13.05 ±0.02 


Si IV 


16.2±4.0 


12.55 ±0.08 


Fell 


(0, 13.9) 


13.31 ±0.15 


SII 


13.9(fixed) 


< 13.37 


CII 


12.7± 1.5 


13.91 ±0.03 


CIV 


13.9(fixed) 


< 12.55 


NI 


13.9(fixed) 


12.79±0.12 


N V 


13.9(fixed) 


< 12.53 


O VI 


(6.9, 14.5) 


13.01 ±0.16 


OI 


10.7±3.8 


13.85 ±0.05 


HI 


13.9(fixed) 


16.84^ 



NOTE. — b is Doppler width of the absorption line, and N is the ionic 
column density. Uncertainties are statistical except for Nm, which includes 
wavelength calibration uncertainties of FUSE observations. See text for de- 
tails. 

Dopplerwidth: v^ C) =-60kms- 1 ,^ I I vc, =6.82x 10 19 cm" 2 , 

and b ( ™ C) = sflo = 17.5 km s" 1 (Fig. __i)- The IVC can also 
be resolved in lines of S n A1259.529 (COS), Si n A1020.699 
(FUSE LiFl a), and O i A929.517 (FUSE SiC2a) - see Fig- 
ures. |2] and [3] We shift the ultraviolet spectra to align their 
IVC velocities with the 21 cm data. 

In the shifted spectra, an HVC at -131 km s" 1 is visible in 
Si ii, O i, N i, and C n at > 3er significance, and in Si iv, 
O vi, and Fe n at 2-3 a significance. All line centroids are 
well aligned within ~ 15 km s" 1 (Figs. [2] and [3]). The HVC, 
however, is not visible in the 21 cm data (Fig.|2^). 

3. HVC MEASUREMENT AND RESULTS 

We use the model developed in l Yao et ail J 20 1 0) to measure 
absorption of the HVC. This model uses ionic column density 
(A/) and Doppler width (bo) as fitting parameters and adopts 
a Voigt function to approximate the absorption-line profile. 
This model can be used to fit individual absorption lines or 
to jointly analyze multiple absorption lines. We include in- 
strumental line spread functions (LSFs) in our spectral fit by 
using the wavelength-dependent LSF of COS observations^ 
and a Gaussian with FWHM 20 km s" 1 across the entire wave- 
length range of FUSE observations. Because Si n has multiple 
transitions (Fig. |2), we jointly analyze these lines to reduce 
statistical uncertainty. This method is equivalent to a curve- 
of-growth analysis but propagates errors automatically. The 
consistent Doppler widths measured in Si n and O i (11-14 
km s" 1 ) indicate that non-thermal broadening dominates the 
line width. To measure column densities or upper limits for 
ions whose HVC components are unresolved or whose line 
widths are unconstrained, we fix their line centroids and their 
widths to those of Si n. Our results are presented in Table Q] 
Because Si m and C m in the HVC are severely blended with 
the IVC components, we have not measured their column den- 
sities. 

We are able to measure A/hi in the high Lyman lines of the 
HVC by subtracting out H i absorption from the IVC and low- 
velocity gas. Both components are saturated, but only the 
IVC (-60 km s -1 ) contributes to the absorption in higher Ly- 
man lines at the (-131 km s" 1 ) velocity of the HVC (Fig. [3]). 

(IVC) , (IVC) 



Since we have well-determined values of v HI v< ~\ b 1 ^^', and 
N ( ™ C> from 21 cm data (§ 13, we fix its absorption in the 



2 http://archive.stsci.edu/ 



3 http://www.stsci.edu/hst/cos/performance/spectral_resolution/ 



High-Metallicity HVC along Mrk 421 Sight Line 



3 



^ 1-5 

S 1.0 

H D 0.5 

0.0 
1.5 

x 

2, i- 

" 0.5 
-a 

.8 ?:g 
H.o 

S 0.5 
z 0.0 



I / 



(a) 



(c) 



H I 21 cm 



-200 -100 



si ii 1260 a (f) 



Si III 1206 A :• (i) 



Si IV 1393 A 




I 1302 A 




VI 1031 A 




100 



-200 



-100 

v,. 



100 

(km s" 1 ) 



-200 -100 



100 



FIG. 2. — LAB 21-cm emission data and UV absorption lines observed along Mrk 421 sight line. Thick black curves are the spectral fit, and blue curves indicate 
local "continuum" of HVC absorption. Vertical dotted and dashed lines indicate centroids of IVC and HVC at —60 km s and —131 km s _1 , respectively. See 
text for details. 
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FIG. 3. — Same as Fig. [2] In panels g-i, data marked with red are spectral ranges used in constraining the HVC H I absorption, and the solid blue, dotted, dashed, 
and dot-dashed curves indicate absorption of the four emission components at -60, -28, —8, and +18 km s _1 revealed in the Green Bank data (Savage et al. 2005). 
Lyr? to LyA are used in our spectral analysis, but only three Lyman lines are plotted for demonstration. See text for details. 



total H i spectrum and fit the remaining absorption as the 
HVC contribution. The deuterium (D i) absorption lines lie 
-82 km s" 1 blueward of the corresponding H i lines. Thus, 
D i in the IVC lies at -142 km s" 1 and gives a small amount 
of contamination to the HVC H i absorption at -131 km s -1 . 
We model the D i absorption by adopting the abundance in 
the local ISM, (D/H) = 1.6 x 10" 5 dLinskv et al.lll993l) and 
assuming that thermal broadening dominates its line width, 
^(ivc) _ ^vc) i^2. To minimize the uncertainties caused by 
line saturation, we only use blue-rising wings of five high- 
order transitions of H i, from Lyij to LyA (Fig. [3). We jointly 
fit these five H i lines to reduce statistical uncertainties, allow- 
ing their positions to vary by ± 10 km s" 1 around the expected 
HVC velocity to account for possible systematic wavelength 
drift. In the joint fit, we set b[f, VC) = C) . The resulting col- 
umn density, logA^ HVC) = 16.84!° f 3 (TableQ) is much better 
determined than the upper limit found from the 2 1 cm data, 
logA^ HVC) < 18.38 (3a), integrated over a velocity range of 
23 km s" 1 , the FWHM of the Si n absorption line (TableQ). 

4. DISCUSSION AND SUMMARY 

The Mrk 42 1 sight line probes an interesting region of HVC 
Complex M. It lies approximately 4° away from the MI, Mil, 
and Mill clouds a nd passes through the IVC IV26 (Fig. [Q 
see a lso Fig. 1 in ISavage etaD 120051 and Fig. 5 in [Wakker 
120011) . An apparent (LSR) velocity change from -112 km s" 1 
at / = 170° to -117 km s -1 at I = 160° was observed in pre- 



vious studies of MI (Wakker & van Woerderifl991l) . The de- 
tected velocity (-130 kms" 1 ) along the Mrk 421 sight line 
is inconsistent with this trend, but it fits into a larger scale 
( > 20°) picture in which the velocities of MI increase toward 
higher longitudes and lower latitudes, connecting with Mill 
whose velocities increase toward lo wer longitudes and higher 
latitudes (Fig.Q) also see Fig. 5 in Wakker 2001). Thus, the 
sight line likely probes the high velocity end of Complex M. 

The HVC along the Mrk 421 sight line has been sug- 
gested in H i, C ii, and C m absorption lines in the FUSE 
spectrum. Savage et al. (2005) measured a surprisingly large 
value, Nl 



(HVC) 
HI 



4.8!^ x 10 1 * cm~ z (logAhi = 18.68 ±0.38) 
in a component fit to the saturated H i Lyman lines Ly/3-LyS, 
nearly 70 times higher than that (Table Q) obtained in this 
work. To examine possible causes of this discrepancy, we 
superpose the amount of the four veloc ity components re- 
vealed from the Green Bank observation (ISavage et al.ll2005l) 
on the H i absorption spectra (Fig. |3g-i). Clearly, only the 
blue wing of the IVC at -60 km s" 1 affects the A^ HVC) deter- 
mination. To further examine possible effects of beam smear- 
ing, we extracted the LAB emission data within 1° around 
the Mrk 421 sight line. Unlike the other three emission com- 
ponents, whose brightness temperatures can vary by a fac- 
tor of five, Nl 



(IVC) 



increases smoothly toward large longitudes 

C C) )^!0% 

among the four beams closest to the Mrk 421 sight line. Tak- 
ing this variation into account in our analysis, we find that 



and high latitudes, with fluctuations AN^ ^ / (N^ 
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the change in A^ VC) is negligible compared to that caused 
by the systematic uncertainty of the wavelength calibration. 
We conclude that our result is not affected by beam smearing, 
although we cannot completely rule out variations of the IVC 
absorption on angular scales < 10 arcmin, which could signif- 
icantly change our results. We speculate that the apparently 
large A^ VC) obtained by Savage et al. (2005) is caused by the 
severely underestimated m their fit. 

The measured A^ VC) and A^ VC) give a metallicity in the 
HVC along the Mrk 421 sight line of [0/H]= 0.32^ or 
0.85-3.5 Zq. The large uncertainty is caused by systematic 
uncertainty in the wavelength calibration introduced in our 
joint analysis of the H i Lyman series (§|3). Because O i and 
H i have nearly the same ionization potential and are coupled 
by resonant charge exchange, this measurement accurately 
reflects the gas-phase oxygen abundance and is not subject 
to ionization correction. In measuring A^ vc \ we assumed 

C C) = C °/V2 and *C C) = *C Q and ad °P ted the local " 
ISM deuterium abundance, (D/ H) = 1.6 x 10~ 5 , o ne of the 
lowest interstellar abundances dRobert et al.l I2000b . Adopt- 
ing higher values of for D/H or b D i will result in a smaller 
iVfH . Therefore the inferred oxygen abundance should be 
regarded as a conservative lower limit. Ionization correc- 
tions need to be applied in calculating the abundances of Si, 
S, or Fe. For a photoionized plasma with hydrogen num- 
ber density «h = 0.1 cm" 3 and logA^i < 17.5, the (logarith- 
mic) ionization corrections are 1.30, 1.42, and 0.90 for Si n, 
S n, and Fe n, respectively (Mark Giroux, private communi- 
cation, 2010). Taking the best-fit values (Table[T|l, we obtain 
[X/H]= -0.6, < 0.22, and 0.07 for Si, S, and Fe. Owing to the 
unconstrained «h, the systematic uncertainty could be high 
in the inferred [Si/H], [S/H], and [Fe/H]. However, this ex- 
ercise results in consistently high metal abundances, except 
for a slightly lower value for Si. Within uncertainties, the 
high metallicit y is consistent with v alues 0.4—1.8 Zp, in MI 
dWakkerll200il) . > 0.5 Z Q in Mill jRvans et alj[l99l. an d 
the near-solar metallicity of IV3 (Spitzer & Fitzpatrick 1993). 
The consistency also suggests that the HVC toward Mrk 421 
is associated with Complex M and the lower velocity IV Arch. 



The measured oxygen abundance provides a new lower limit 
to the metallicity of Complex M. 

The high metallicity, together with the upper limit on dis- 
tance of other parts of Complex M, shed light on its ori- 
gin. Unlike the highly ionized HVCs, in which the high ions 
(O vi and C iv) are significant or even dominant, the HVC 
toward Mrk 421 is dominated by low ions (O i, C n, S n; 
Table [TJ. This suggests distinct ionization conditions of this 
HVC comp ared to the highly ion i zed, low H i column den- 
sity HVCs (ISembach et all 120031: ICollins et al.l l200l l2005t 
IFox et al J 120061: iLehner & Howkl 120 1 01) . Such mildly ion- 
ized, low H i column density HVCs have recently been sur- 
veyed by several groups ( e.g., Ric hter et al.ll2009l IShull et al.l 
120091: ICollins et alJUool . Richter et al. (2009) also pro- 
vided ionization modeling for the HVCs with different dis- 
tances. Substituting A^ VC) and A^ vc) (Table Q) into their 
model and assuming no depletion of Si into dust grains, we 
obtain a gas density «h = 3.4 cm" 3 and absorption path length 
L = 0.007-0.02 pc for Complex M. The near-solar (likely su- 
persolar) abundance, short path length, and large gas density 
clearly rule out the WHIM model and disfavor the halo model 
(§03 for the HVC. Complex M likely traces a supernova shell 
that has not been ejected too far above the Galactic plane and 
is now falling back to the Galactic disk. 

In summary, we have presented an interstellar high veloc- 
ity cloud detected in ultraviolet spectroscopic observations of 
Mrk 421, which represents the highest velocity end of the 
HVC Complex M. The combination of HST-COS and FUSE 
observations enables us to obtain a lower limit, > 0.85 Zq, 
to the metallicity of the Complex. The high metallicity and 
the short distance are consistent with Complex M being the 
returning gas of a Galactic fountain. 
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